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Abstract. We study a class of extended quintessence cosmologies where the scalar 
field playing the role of the dark energy is exponentially coupled to the Ricci scalar. 
We find that the dynamics induced by the effective gravitational potential in the 
Klein-Gordon equation dominates the motion of the field in the early universe. The 
resulting "i?— boost" trajectory is characterized by a kinetic dark energy density, given 
by [Spmnroi^ + z)] [32pr0 ^ jBD o] ~ , where lojbdo, pro and pmnro are calculated at 
present, and represent the Jordan Brans Dicke parameter, the density of relativistic 
matter and of those species which are non- relativistic at redshift z, respectively. 
We show that such a trajectory represents an attractor, equivalent to a tracking 
solution with equation of state w ~ —1/3, providing a large basin of attraction for 
the initial dark energy density regardless of the properties of the potential energy 
yielding acceleration today. We derive the up to date constraints from Big Bang 
Nucleosynthesis (BBN) on the present scenario, and we show that they are largely 
satisfied for interesting trajectories of the dark energy scalar field in the early universe. 
We compute the cosmological perturbation spectra in these cosmologies. For a fixed 
value of ujjBDo, the projection and Integrated Sachs- Wolfe effects on the cosmic 
microwave background anisotropy are considerably larger in the exponential case with 
respect to a quadratic non-minimal coupling, reflecting the fact that the effective 
gravitational constant depends exponentially on the dynamics of the field. 
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1. Introduction 

According to tlie present cosmological observations, the universe today is nearly 
geometrically flat, with an expansion rate of about 72 km/sec/Mpc, and structures 
grown out of a primordial linear spectrum of nearly Gaussian and scale invariant energy 
density perturbations. About 5% of the critical energy density is made of baryons, 
while 20% is thought to be composed by particles interacting at most weakly with them 
(Cold Dark Matter, CDM). The remaining 75% is some sort of vacuum component, 
the dark energy, with negative pressure acting as a repulsive gravitational force, and 
responsible for a late time cosmic acceleration era. Several independent observables 
support the picture above, i.e. type la supernovae (hereafter SNIa P3I2]), the Cosmic 
Microwave Background (CMB) anisotropies (see jS] and references therein), the Large 
Scale Structure (LSS, see jlj and references therein) and the Hubble Space Telescope 
(HST 5J). 

While a dark matter component is generically predicted in the super-symmetric 
extensions of the standard model for particle physics, the dark energy component poses 
severe questions to the whole picture (see jHl [Z| and references therein). Specifically, a 
reliable dark energy candidate should explain why the present amount is so small with 
respect to the fundamental scale (a fine tuning of 123 orders of magnitude if compared to 
the Planck scale), and why it is comparable with the critical density today (coincidence). 
The Cosmological Constant is affected by both these problems. As an alternative, the 
concept of vacuum energy in cosmology has been generalized in terms of a scalar field, 
named quintessence, which for a broad class of potentials is able to converge to the 
present regime starting from a wide set of initial conditions in the past, thus alleviating 
the fine tuning (scaling or tracking trajectories, 9j). The lack of a clear prediction 
for the dark energy scalar field in any fundamental theory motivated the investigation, 
also phenomenological, of its coupling with other cosmological components, such as dark 
matter [TUl [TT] or gravity [12] . 

This work belongs to the second class. The interest in connecting dark energy 
and gravity is related to the original Brans-Dicke idea (see JH] and references therein), 
and on the temptation to explain the evidence for cosmic acceleration entirely in terms 
of a modification of general relativity. Although such attempts were unsuccessful so 
far, several interesting effects have been discovered in these scenarios, generally called 
extended quintessence models, concerning background evolution as well as cosmological 
perturbations [12], [I1-[2H]- 

A major achievement concerns the early universe dynamics of the quintessence field 
in these scenarios. It was shown that the species which are non-relativistic in the 
radiation dominated era make the Ricci scalar diverging, activating an effective potential 
in the Klein-Gordon equation coming from the non-minimal coupling only (i?-boost, 
|19j). This effect is acquiring a crucial importance recently: as we stressed above, a 
broad class of quintessence potentials allows for a large basin of attraction in the early 
universe, alleviating the corresponding fine tuning problem. On the other hand, the 
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latter condition holds only if the present dark energy equation of state w is sufficiently 
far from —1. As the observations are constraining w to be close to —1 (see e.g. j3]) 
the basin of attraction in the early universe shrinks, threatening the basis of the whole 
quintessence picture [21] • It was recently argued that in these conditions the coupling 
of the dark energy with other entities may be relevant at early times, possibly saving 
the existence of attractors for the initial field dynamics [T^ . 

Despite the broad interest on this subject, most of the analysis so far was carried 
out considering a quadratic coupling between the dark energy scalar field and the Ricci 
scalar. In this paper we work out the first generalization of those claims, providing 
a comprehensive analysis of cosmologies involving an exponential coupling. The latter 
case may be relevant for string cosmology, where the dilaton appears in an exponential in 
the fundamental lagrangian [201, although in the usual formulation the coupling involves 
all terms, not just the Ricci scalar. 

This paper is organized as follows. In Section II we describe the framework of 
exponential couplings in extended quintessence cosmology. In Section III we discuss the 
observational constraints on the coupling magnitude from solar system experiments and 
Big Bang Nucleosynthesis. In Section IV we study the -R-boost trajectories of the dark 
energy scalar field. In Section V we work out the most important effects on cosmological 
perturbations. Finally, in Section VI we draw our conclusions. 

2. Exponential scalar-tensor couplings in cosmology 

A wide class of cosmological models can be described by the action 



where g is the determinant of the background metric, R is the Ricci scalar, (j) is 
a scalar field whose kinetic energy and potential are specified by uj{(f)) and V{(j)), 
respectively. /^Auid includes contributions from all components different from and 
K = SttG^k represents the bare gravitational constant [201 • The usual gravity term 
R/lGnG has been generalized by the function f{(j),R)/2K,, where f{(f),R) may describe 
a coupling between the quintessence field </> and gravity, as well as a pure geometrical 
modification of general relativity, featuring a non-linear dependence on R. Note also 
that the gravitational constant appearing in the Lagrangian and the one measured in 
Cavendish like experiments differ by corrections being negligible in the limit ujjbd ^ 1 
PU] . The Lagrangian in has been considered and analysed in its full generality 
in a cosmological context including linear perturbations [3T]; this framework has been 
exploited in several works [121 dZl UHl 1211 1211 l2Ej , leading to the formulation of the weak 
lensing theory [2Zj- 

In this work we consider a simple class of extended quintessence models with 
exponential coupling 




(1) 




(2) 
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where ^ is a dimensionless coupling, 0o is the present value for the field and mp 
the Planck mass irip = l/\fG in natural units. Note that the constant 0o has been 
introduced to make explicit that at present -F(0o) = I/SttG. 
The Jordan-Brans-Dicke parameter in this scenario is 



'dF\ ^ Stt 
UJJBD = \ — \ =7:jexp 



mp 



(3) 
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Note that at present ujjbdo = Svr/^^. 

We give now a brief overview of the relevant equations describing the cosmological 
expansion assuming a Friedmann Robertson Walker (FRW) background metric, since 
it is relevant in the following. In Section V we describe the effects on the cosmological 
perturbations; however, we do do not write down explicitely the perturbation equations 
here, which can be found elsewhere (see e.g. 1^1211) written for a generic function 
F. The most relevant differences with respect to ordinary cosmologies are represented 
by the time variation of F, which plays the role of a varying gravitational constant in 
the Poisson equation. 

For fiat cosmologies, the line element can be written as ds"^ = a^(r)(— dr"^ + 5ijdx^dx^) , 
where a(r) is the scale factor, r represents the conformal time variable, related to the 
cosmic time t by the transformation dt = a(r) dr. 

The expansion and field dynamics are determined by the Friedmann and Klein Gordon 
equations 

= (^) ' = ^ {a'pfim, + l^' + a'V- mP) , (4) 

-^ + 2n^=-—R-a'V^, (5) 

where dot means derivative with respect to the conformal time r and is the derivative 
of the quintessence potential with respect to 0. For our analysis, it is relevant to write 
down the explicit form of the Ricci scalar in terms of the cosmological content 



-4 



F UF" 



-pfiuid + ^Pfiuid + ^-4V + 3 — + 2 



(6) 



where p fluid and p fluid are the energy density and pressure that need to be summed up 
over all possible cosmological components but 0. Note also that, as it was shown in |31j . 
the usual conservation equations pi = — 37i(pj + Pi), where pi and Pi are respectively 
the energy density and pressure of the i-th component, hold for all species but 0. 



3. Observational constraints 



In this Section we briefiy review the constraints on scalar-tensor cosmology coming from 
the solar system physics [22] and derive the up to date bounds from the BBN. 
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The time variation of the gravitational constant is related to the coupling function F 
by the expression 

dG _ IdF 

which does not take into account corrections due to the effective gravitational constant 
in Cavendish like experiments which are small in the limit ujjbd ^ 1 [20] ■ Both local 
laboratory and solar system experiments constrain the ratio in (|7j) to be less than 
10~^^ per year jHS]. An independent bound is given by the effects induced on photon 
trajectories by a varying gravitational constant, which was recently greatly strengthened 
by the Cassini solar system probe [33j: at two sigma, the new constraint is 

Cc-JBDO >4xl0^ (8) 

It may not be straightforward to extend these limits, obtained on solar system scales, 
to cosmology, as it was pointed out recently [31]. Indeed we probe regions well within 
our Galaxy, which is a self-gravitating virialized system: in physical theories where 
fundamental constants vary, the latter may acquire local values which are different 
from their large scale, cosmologically effective ones. For this reason, cosmology is 
likely to become a source of constraints for the underlying theory of gravity ^35^, in 
a complementary way with respect to the solar system. 

The BBN is indeed the only direct cosmological probe for the field value in the 
radiation dominated era. As it is well known, the amount of light nuclides produced 
when the photon temperature was in the range 0.01 -^ 10 MeV is rather sensitive to 
the value of the Hubble parameter during that epoch, as well as to its time dependence. 
In particular the ^He mass fraction Yp strongly depends on the freeze-out temperature 
of weak processes which keep neutrons and protons in chemical equilibrium. Changing 
the value of Ti affects the neutron to proton number density ratio at the onset of the 
BBN, which is the key parameter entering the final value of Yp and more weakly in 
the Deuterium abundance. For a fixed baryon density parameter ujb = VLi^h? in the 
range suggested by results of the Wilkinson Microwave Anisotropy Probe (WMAP), 
ujb = 0.023±0.001 inni both nuclei yields are in fact monotonically increasing functions 
of H. For recent reviews on BBN see |37 ] l38 | l39]. 

As we read from ^ the introduction of the field changes the value of 7i in two 
different ways. First of all there is a shift in the effective gravitational constant by a 
factor exp[— ^(0 — (f)o)/'mp)]. Second, it gives an extra contribution to H. By using the 
coupling defined in (0), we can solve ^ with respect to H, getting 



n 



\ {p fluid + P<i>) exp -J ) + J— 

\ V V ^JBDo rnp J AujjBDo V 4cu. 



JBDO 



.(9) 



In the following we shall assume that the Universe is indeed radiation dominated during 
the BBN epoch, so that p<^ can be neglected with respect to p fluid which receives 
contributions from photons, neutrinos and e^. Moreover the last two terms in Equation 
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Q are suppressed by 1/ y/uJjBDo- In this case the squared Hubble parameter takes the 
simple form 



H = — = — -exp 



P fluid ■ (10) 

6 [ TJip 

In this limit we see that by comparing the theoretical values for the light nuclei 
abundances with the corresponding experimental determinations we can constrain the 
value of the effective gravitational constant during the BBN epoch, i.e. the quantity 
^(0 - (t)o)/mp. 

Before discussing this issue let us briefly summarize the present knowledge of the 
BBN both from the experimental and theoretical points of view, see e.g. j3H| for a review. 
The most accurate measurement of primordial Deuterium number density, normalized 
to Hydrogen, X/j, is obtained from DI/HI column ratio in QSO absorption systems at 
high redshifts. The most recent estimate jlHI gives 

X^^^=(2.78^°:^^)-10-^ , (11) 

and is in good agreement with the theoretical expectation for a standard scenario with 
three active neutrinos and a baryon density given by the WMAP result the latter 
is 

X'^ = (2.44+0:i?)-10-5 , (12) 

where the (la) theoretical uncertainty accounts for both the propagated error due to the 
several rates entering the BBN nuclear reaction network as well as the 5% uncertainty 
on ujh- On the other hand the ^He mass fraction Yp obtained by extrapolating to zero 
the metallicity measurements performed in dwarf irregular and blue compact galaxies 
is still controversial and possibly affected by large systematics. There are two different 
determinations ^TJ |12] which are only compatible by invoking the large systematic 
uncertainty quoted in [Hj 

Fp = 0.238±(0.002),i„t±(0.005),j,, , (13) 
Yp = 0.2421±(0.0021),t,t . (14) 

Both results are significantly lower than the theoretical estimate jHEj 

Fp*^ = 0.2481±0.0004 , (15) 

where again we use cj;, = 0.023 and the small error is due to the uncertainty on the 
baryon density and, to a less extent, to the error on experimental determination of the 
neutron lifetime. As in [3Hj we use in our analysis a more conservative estimate for the 
experimental "^He abundance obtained by using the results of [43j 

Y^'^'P = 0.245±0.007 . (16) 

We finally recall that the Spite plateau value for ^Li found in PopII dwarf halo stars is 
smaller by a factor 2^3 with respect to the standard BBN prediction p^l39j. In view 
of this problem, suggesting a large depletion mechanism of primordial ^Li, we do not 
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consider this nuclide in the following. 

To bound the value of the quintessence field we construct the likelihood function 
with 

x'm-<Po)/mp] = 



[17) 



^exp 



w,.. 



xf[e(0-0o)/mp] 



[18) 



E [Xf[e(0-0o)/mp]-xr 

i,i=D,4He 

The proportionality constant can be obtained by requiring normalization to unity of the 
integral of C, and with Wij we denote the inverse covariance matrix 



+ '^i,exp^ij] 



(19) 



where o"j is the uncertainty in the experimental determination of nuclide abundance 
Xi and afj the theoretical error matrix. We also consider the two likelihood functions 
for each of the two nuclei to show how at present the Deuterium and ^He can separately 
constrain the value of the effective gravitational constant at the BBN epoch. 

It is important to stress that we consider the simplest case of a constant value 
of (p during the whole BBN phase. This is justified a posteriori by the discussion in 
the next Section; although the dynamics of the field in the radiation dominated era is 
cosmologically relevant, it is too small to provide any change of during BBN which is 
significant compared with the bounds we derive here. It is also worth mentioning that 
in view of the possible systematics affecting mainly the experimental determination of 
'^He it is unfortunately impossible at present to use BBN to get detailed constraints on 
the time evolution of the quintessence field during this phase, while it can only provide a 
conservative bound on the largest (or smaller) values attainable by C,{4> — <f)o)/'mp during 
the nuclei formation era. 




-0.6 -0.4 -0.2 0.2 



Figure 1. The behavior of the lilcelihood functions versus the parameter ^{(j)—(j)o)/mp 
for LJb — 0.023. The three curves refer to the D (dashed curve), '^He (sohd) and the 
combined D-|-^He (dotted) analysis discussed in the text. 



Our results are summarized in Figure ^ where we show the likelihood contours 
obtained using Deuterium, "^He and finally, combining the two nuclei. As expected the 
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^He mass fraction gives the most stringent constraint since it is much more sensitive to 
the value of the Hubble parameter with respect to Deuterium. In particular using ^He 
only we get 

- 0.13 < ^(</) - (/)o)/mp < 0.20 , 95% C.L. (20) 
while adding the information on D abundance does not change this result significantly: 

- 0.14 <^(0-0o)/mp< 0.17 , 95% C.L. . (21) 

We see how a clear understanding of the role of systematics in the Yp measurements 
would have a large impact on further constraints on the value of the effective 
gravitational constant, at least in the minimal BBN scenario we are here considering, 
with no other extra parameters, such as extra relativistic degrees of freedom. 

We close this Section by recalling that for a constant value of (p during the BBN 
epoch, its effect in renormalizing the gravitational constant can be also conveniently 
recasted in terms of the effective number of relativistic species contributing to the total 
energy density, as it was noticed in earlier works fTfl I44j . We parameterize the value 
of the Hubble parameter at the onset of BBN and before the annihilation phase in 
terms of an excess (or deficit) of relativistic degrees of freedom AA^ 

where p-y is the photon energy density and = 10.75 the standard value obtained by 
summing over photons, and neutrinos. By using (fTO|l it is straightforward to get the 
simple relation 

5i^ = _,og(l + i^Aiv) . (23) 

mp V 43 / 

This relation can be used to recast the bounds on the effective number of neutrinos, 
which are routinely used in the literature, in terms of a constraint on C,{4> — (f)o)/mp 
for the exponential coupling considered in this paper. As we show in the next Section, 
interesting values of the field at Nucleosynthesis largely satisfy the bounds represented 
in Figure n 

4. i?-boost 

In this Section we derive the dark energy dynamics before the onset of acceleration. The 
central point is that extended quintessence scenarios may possess attractor trajectories 
in the early universe, generated only by the non-minimal coupling, as it was pointed 
out recently jT2]; that allows to remove the fine tuning in the early universe, even if the 
potential is constrained to be almost flat in order to be consistent with the observations; 
the latter occurrence indeed limits the capability of minimally coupled quintessence 
models to avoid the fine tuning problem [221 • The attractors in minimally coupled 
quintessence models were called scaling jSj or tracking P solutions; the corresponding 
trajectories in extended quintessence coming from the non-minimal coupling are called 
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i?-boost solutions 19 . In the following we derive the i?— boost for exponential 
coupling. We shall write an analytic expression for that, both in matter dominated and 
radiation dominated eras, which is manifestly independent on the initial conditions; 
that represents an important new aspect with respect to earlier works in which 
the initial field value appeared explicitely in the i?— boost energy density. Moreover, 
exploiting a perturbative analysis [H], we demonstrate that such solution represents an 
attractor. Also, we shall see that the energy density corresponding to the i?— boost 
depends on the value of ujjbd- 

We shall compare what we find here with previous models of extended quintessence, 
based on a quadratic coupling between the dark energy and the Ricci scalar, and with 
a case of minimally coupled quintessence (QCDM). In order to satisfy the constraint 
(jHl) the constant ^ is obtained through by fixing uijbdo = 10^. The remaining 
cosmological parameters are chosen consistently with the concordance model (see e.g. 
PE]). The present dark energy density is 73% of the critical density, fl^ = 0.73, 
with a Cold Dark Matter contribution of ^Icdm = 0.226, three families of massless 
neutrinos, baryon content = 0.044 and Hubble constant Hq = 72 Km/sec/Mpc. The 
quintessence potential is an inverse power law 

y(<^) = ^ ' (24) 

with a close to zero in order to have close to —1 at present (see e.g. and 
references therein). We choose the initial conditions at 1 + z = 10^; in typical runs, with 
= —0.9 at present, the value of the field today is 0o = 0.35mp; the dynamics induced 
by the potential V and the i?— boost makes <^ 0o at early stages; the condition on 
^jBDo sets the value of the coupling constant ^ = 1.6x10"^. As a consequence, the 
BBN bounds (pUjl or (j2H) are largely satisfied. 

4.1. Radiation dominated era (RDE) 

The crucial point is that in (jH)) the term p fluid— 'ip fluid gets no contribution from radiation 
but, as long as one or more sub-dominant non-relativistic species exists, the Ricci scalar 
gets a non-zero contribution that diverges as 1/a^, thus becoming increasingly relevant 
at early times 

1 PmnrO 

The parameter Pmnro above indicates the energy density of the species which are non- 
relativistic at the time in which the dynamics is considered. As a consequence, the term 
proportional to R in the Klein Gordon equation acts as a new source on the quintessence 
field besides the true potential V{(f)); the resulting motion is the i?— boost. 
With our choice of the coupling Q the Klein Gordon equation ^ for the field can be 
written as 

4> + 2n^-l^^^ + a% = , (26) 
2 mp a 



Rr^±[^ fora^O . (25) 

77' 
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where we have used (0) to write the Ricci scalar in terms of the total energy density, 
neglecting the remaining terms in © since they get negligible for a sufficiently small. 
For our purposes here we neglect the departure from general relativity in the RDE, i. e. 




Figure 2. Energy density of matter (solid), radiation (heavy dashed) and 
for minimal (dot dashed), quadratic (heavy dotted) and exponential (heavy solid) 
coupling. Also shown the analytical behavior of the R-boost during radiation 
dominated era (dashed) and matter dominated era (dotted). 



we assume F ~ I/SttG and we ignore all the contributions in the Friedmann equation 
except radiation. The above assumption is justified because we aim at deriving the 
-R— boost dynamics at first order in 1/uijbd'- a general analysis would yield corrections 
of higher order in 1/ujjbd in the effective gravitational potential appearing in the Klein 
Gordon equation. The expansion is given by 

a ~ , (27) 



where C = ySTrGpro/S and p^o represents the cosmological radiation density calculated 
at present; since 7i ~ 1/r, the equation of motion for (j) reduces to 

0+-0-^^^=O , (28) 
r Imp Ct 

which is solved by 

c I Q 

-{r - Tbeg) + (pbeg , (29) 

where (pbeg is the initial condition for 0. This solution corresponds to the slow roll phase 
which starts when the cosmological friction and effective gravitational potential effects 
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in the equation of motion are comparable, yielding the i?— boost equation 

• a^RdF , , 

2W.— - . (30) 

In this phase, the energy of the quintessence field is dominated by its kinetic contribution 

3 pmnrO 1 




-(1 + ^r • (31) 

from which it is possible to derive that the R-boost solution is equivalent to a tracking 
one with equation of state —1/3. It is interesting to note that at first order inl/ujBD the 
i?— boost energy density is related only to the present value of the Jordan-Brans-Dicke 
parameter. In Figure |21 we have plotted the energy density of the various cosmological 
components: besides matter and radiation we can observe the behavior of (j) for minimal 
and extended quintessence, starting from zero initial kinetic energy. In the minimal 
coupling case the field behaves nearly as a cosmological constant until the true potential 
V starts to be relevant. In the extended quintessence case, instead, the field accelerates 
and soon enters the i?— boost phase. In Figure 01 we show the absolute values of the four 
terms in the Klein Gordon equation (f^H|) . The potential term starts to be dominant only 
for z < 10^; the friction term is zero at the beginning and then increases, joining the 
i?— boost. The quintessence field accelerates until the sign inversion in occurs; then 
the field accelerates again for z ~ 10^ when the potential V starts to have a relevant 
effect. The timing of the different phases of the trajectory depends on the details we 
have fixed, but the general behavior is model independent. 
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Figure 3. Absolute values of the four terms in the Klein Gordon equation H28|) for an 
exponential coupling. 
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4-2. Matter dominated era (MDE) 

In the matter dominated era, equation ()30|1 still holds but the expansion parameter has 

a different behavior with time 
2 

a = -nGpmnroT ■ (32) 

As in the RDE case, we neglect the departure from general relativity to get the i?— boost 
dynamics at the lowest order in I/ojjbd- The i?— boost equation is 

1 3 e 



IGtt r 

which is solved by 



mp , (33) 



3 r 

= T^i^P log ^ (t>beg , (34) 

167r 

also shown in Figure |21 As a consequence, the behavior of the kinetic energy of the field 
changes, too 

2 . ^ 

-Pmnr0(l + zf . (35) 




32 UJjBDQ 

If we now look at Figure El again, we notice that the i?— boost has a bigger effect on p in 
the exponential case with respect to the quadratic coupling, also shown in the figure. In 
both cases receives a strong kick, which determines a major change in the dynamics 
with respect to the minimally coupled scenario. However, for an exponential coupling 
the path of the field departs from the standard one earlier with respect to the case 
of a quadratic couphng, reflecting the fact that the exponential coupling enhances the 
departure from general relativity as it depends exponentially on the field dynamics. We 
shall come back to this issue in the next Section. 



4-3. Stability 

Following the method exploited in , we now show that the i?— boost solution above is 
an attractor. We look for scaling solutions of the Klein Gordon equation, i. e. solutions 
in which the energy density of quintessence field scales as a power of the scale factor 



\dt 

in the RDE and 



oc a , (36) 



{-ft) • P^' 

in the MDE, as found in (j31l The behavior of the background energy density is 
also given by a power of the scale factor. In particular 

Pr OC , pm oc , (38) 

for radiation and non-relativistic matter respectively. As far as a{t) is concerned, 
neglecting again the corrections of the order I/ojjbd induced by 1/F in front of the 
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right hand side of the Friedmann equation (@)) as well as the other dark energy terms, 
one has 



^2/3 . (40) 



The corresponding time dependence of the field (p is given by 



RDE 



t. 



(41) 



<^e^^^ = 0*log(^) , (42) 

where the subscript stands for a given reference time and the subscript 'e' reminds 
us that this is the R— boost exact solution. 

We verify that the i?— boost solution is an attractor by rewriting the Klein Gordon 
equation with the change of variables given by 

u = ^ ' (43) 

(Pe 

and 

r = log(i) . (44) 

During RDE, we obtain 

2u" + 3m' + M - 1 = , (45) 

where the prime stands for the derivative with respect to f. As we immediately see, 
this equation admits a critical trajectory for u = 1 and u' = 0, corresponding to the 
i?— boost. Linearizing equation pSjl by choosing small exponential perturbations around 
the critical point [u = 1 + e^'^) and solving for the eigenvalues Ai^2 we get 

Ai = -1 ; A2 = . (46) 

The fact that the eigenvalues are real and negative shows that the generic perturbation 
will be suppressed with time, thus flattening the trajectory on the critical one = (pe- 
Similarly, during the MDE, the Klein Gordon equation (jSJ can be rewritten as 

u" + u = , (47) 

which again admits m = 1 as a solution. In this case, by substituting u = u + e^'^ in the 
equation and solving for the A eigenvalues, we find the two values Ai = and A2 = — 1, 
independently of the value of u. As a consequence, the generic solution of the Klein 
Gordon equation will tend to m = constant. The constant is fixed to be equal to 1 by 
the initial conditions for the MDE, corresponding to the final regime in the RDE, where 
the trajectory joined the i?— boost solution represented by (pe- 

Summarizing, the simple choice of an exponential as a coupling between dark energy 
and Ricci scalar in the Lagrangian leads to the existence of an attractive i?— boost 
solution in both the matter and radiation dominated eras. 
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Figure 4. The behavior of if as a function of the redshift in three cosmological 
models, extended quintessence with exponential coupling (solid), quadratic coupling 
(dashed) and QCDM (dotted). 

5. Effects on cosmological perturbations 

In this section we analyze the effect that the exponential coupling in (Q) has on the 
spectra of cosmological perturbations. We find that these effects are qualitatively 
consistent with those corresponding to a quadratic coupling, but quantitatively different. 
For our purposes, an analysis based on the CMB power spectrum only is sufficient. 
The scalar perturbations are Gaussian and described by a scalar power spectrum with 
spectral index n = 0.96, and no tensors, consistently with the cosmological concordance 
model, see jSH] and references therein. 

Since the constant ^ is chosen to be positive and the field is smaller in the past 
with respect to (po, F < l/SvrG in the past. As a consequence, our model describes 
a cosmology in which the gravitational constant (and thus the Hubble parameter) is 
higher in the past than in the QCDM case. It follows that for a fixed ujjbdq-, we expect 
a larger amplitude of the effects induced by the behavior of F for the exponential case 
with respect to the quadratic coupling: indeed, the exponential is sensitive to the field 
dynamics also at the linear level, which dominates for small values of 0/0o ~ 1- This 
can be seen in Figure |31 where we plot H{z) for three cosmological models: two of them 
represent extended quintessence with the same value of ujjbd at present, but featuring 
an exponential and quadratic coupling, while the third one is the QCDM case. For the 
first two cases, we have set ujjbdo = 50; we stress that we have chosen this small value 
with respect to the existing bounds from solar system in order to highlight the 
differences in the models considered. The value of Hq is the same, but H{z) are quite 
different functions in the three cases, in particular for the exponential case with respect 
to the QCDM cosmology. 

In Figures El and El we plot the total intensity and polarisation power spectra of 
anisotropics for the models considered in Figure ^ Two features are immediately 
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evident. First the different amplitude in tlie tail at low multipoles, and second a 
projection difference in the location of the acoustic peaks. The spectra have been 
normalized to the amplitude of the first peak, which roughly corresponds to fixing 
the amplitude of the signal at last scattering. Both the effects come from the different 
behavior of the coupling F with redshift. The projection feature simply follows from 
the difference in the curves in Figure |3 Indeed, smaller (higher) values of project 
the CMB power spectrum onto smaller (larger) angular scales in the sky. The power 
at low multipoles is modified through the Integrated Sachs- Wolfe (ISW) effect, which 
is sensitive to the change of the cosmic equation of state at low redshift (see jlH] and 
references therein). The larger is that change, the larger is the ISW power. 

Summarizing, the different form of the coupling function F determines a relevant 
difference of amplitude in the effects induced by extended quintessence models on the 
cosmological perturbations. Larger departures from general relativity in the dynamics 
of the field in the early universe correspond to larger effects on the perturbations. The 
latter aspect has to be taken into account when constraining different theories on the 
basis of the present observations. 




Figure 5. CMB angular total intensity power spectra for a QCDM cosmology 
(dotted), quadratic (dashed) and exponential coupling extended quintessence (solid) 
with ujjBD = 50. The spectra are in arbitrary units, normalized to 1 at the first 
acoustic peak in total intensity. 



6. Conclusions 

We made a general analysis of quintessence cosmologies where the dark energy scalar 
field is exponentially coupled with the Ricci scalar. With respect to the case of a non- 
minimal quadratic coupling, we find relevant new results concerning both background 
expansion and perturbations. 

The dynamics of the field in the early universe is dominated by the effective 
potential in the Klein-Gordon equation, entirely due to the non-minimal interaction. 
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Figure 6. CMB angular polarisation power spectra for a QCDM cosmology (dotted), 

quadratic (dashed) and exponential coupling extended quintessence (solid) with 
^JBD = 50. The spectra are in arbitrary units, normalized to 1 at the first acoustic 
peak. 

The presence of non-relativistic species, although sub-dominant with respect to 
radiation, makes the behavior of the Ricci scalar actually diverging in the early universe, 
determining the dynamics of the field. The whole trajectory, the i?-boost, corresponds 
to the slow roll of the field on the effective gravitational potential. 

Interestingly, in the present scenario the i?-boost looses any sensitivity on the initial 
conditions of the field, while in the case of a quadratic coupling the latter enters into 
the expression for the energy density along the trajectory. The i?— boost solution is 
then equivalent to a tracking one with equation of state —1/3, and energy density 
[3pmnro(l + z)f ^2prQUJjBDoT^ , where Pro is the radiation energy density today, pmnro 
the present energy density of the species which are non-relativisitc at redshift z, and 
^.JBDO is the Jordan-Brans-Dicke parameter, also calculated at present. This aspect is 
extremely relevant in dark energy cosmology; indeed, the basin of attraction of the true 
quintessence potential shrinks if the equation of state today gets close to —1 as suggested 
by the observations. In this context, it is important to check whether the coupling of the 
field with other entities, gravity in this case, is able to restore the existence of attractors 
for the initial dynamics of the scalar field. We find that if the coupfing is exponential 
the answer to that question is positive. 

We studied the impact of this scenario on the observations. In particular we derived the 
limits from the Big Bang Nucleosynthesis; the constraints affect both the strength of 
the non-minimal coupling as well as the field value during the nuclide formation epoch. 
For interesting trajectories, those bounds are largely satisfied. 

We have finally considered the Cosmic Microwave Background (CMB) perturbation 
spectra, comparing the effects in the cases of exponential, quadratic and minimal 
couplings. The main effects are due to the time variation of the effective gravitational 
constant; the latter represents an extra-source of dynamics for the Hubble expansion 
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rate. That modifies distances and perturbation growth rate, shifting the angular location 
of the acoustic peaks in the CMB power spectrum, and enhancing its tail at low 
multipoles, by means of the modified dynamics of the gravitational potentials on large 
scales (Integrated Sachs- Wolfe). For a fixed ujjbd at present, a marked difference exists 
between the amplitude of these effects in the two extended quintessence cosmologies 
considered. Indeed, the exponential form enhances the departure from general relativity 
as the field moves, with respect to the case of a quadratic coupling. The reason is just 
the shape of the coupling, for which the motion of the field gets exponentially amplified. 

In conclusion, two main remarks can be made. The first one is the capability of 
generalized theories of gravity to provide attractors for the early universe dynamics 
of the dark energy scalar field, independently of the shape of the potential inducing 
acceleration. The relevance of this issue goes beyond the mere capability of avoiding 
the fine tuning on the initial conditions. Actually, this aspects shows that the relation 
between dark energy and other entities deserves a careful study, providing a rich 
physics even if the present behavior is forced by observations to be close to that of 
a cosmological constant. The second remark concerns the capabihty of cosmology to 
constrain modified theories of gravity. Interestingly the imprint on the cosmological 
perturbations spectra depend on the actual form of the non-minimal interaction. The 
analogy existing between different scalar-tensor theories of gravity cannot be extended 
to their imprint on cosmological perturbations, which has to be done properly case by 
case. This aspect is relevant having in view of the possibility to constrain the theory of 
gravity on cosmological scales and from cosmological observations, in a complementary 
way with respect to what is done via observations in the solar system. 
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